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Abstract 

Gibbs' measures in the Sherrington-Kirkpatrick type models satisfy two asymptotic 
stability properties, the Aizenman-Contucci stochastic stability and the Ghirlanda- 
Guerra identities, which play a fundamental role in our current understanding of these 
models. In this paper we show that one can combine these two properties very naturally 
into one unified stability property. 
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1 Introduction and main results. 

In the Sherrington-Kirkpatrick (SK) model [29J or, more generally, in the mixed p-spin 
model one considers a random process (Hamiltonian) H^(cr) indexed by spin configurations 
cr G = { — 1, + 1}^ given by a linear combination J2 P >i PpHn, p (ct) with some coefficients 
P p > of the independent Gaussian processes, called p-spin Hamiltonians, 

HnM) = N {p-l)/2 9h,...,i p Vii ■ ■ ■ o ip , (1.1) 

l<ii,...,i p <N 

where (^h,...^) are standard Gaussian independent for all p > 1 and all (ii, . . . ,i p ). The 
Gibbs measure Gn corresponding to the Hamiltonian is denned as a random probability 
measure on E^v given by 

G N ((r) = ^ex P H N (cT), (1.2) 

where the normalizing factor Zn is called the partition function. The Gibbs measure Gn m 
fll.2[) is the central object of interest in spin glass models and the answers to many important 
questions follow from the conjectured properties of Gn- These properties can be expressed in 
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terms of various functions of the sample (cr l )i>i from Gn, for example, the normalized Gram 
matrix R = {Ri,v)i,v>i = N' 1 ^ 1 ■ er l )i > i'>i- It is easy to see that knowing R is equivalent 
to knowing Gat up to orthogonal transformations since one can reconstruct Gn from R up 
to orthogonal transformations (as we shall see in the proof of Theorem [1] below) and the 
information encoded in the distribution of R turns out to be sufficient for most purposes 
in the setting of the SK model due to the fact that the Hamiltonian Hn{<t) is a Gaussian 
process and its covariance KHn((t 1 )Hn((t 2 ) is the function of exactly the normalized scalar 
product R12 = N~ l a l ■ cr 2 , called the overlap of cr 1 and cr 2 . Given any limiting distribution 
of the Gram matrix R in the thermodynamic limit iV — > oo, one can use the Dovbysh- 
Sudakov representation ([IQ],[T9]) to define the asymptotic analogue of the Gibbs measure 
as a random probability measure G on the unit ball of the Hilbert space £2 (see [5], [3], 
[15]). This means that in the limit the matrix R is still generated as the Gram matrix of 
the sample from some random measure G. For simplicity, let us assume that this asymptotic 
Gibbs measure G is atomic, 

G = X>*6' (1.3) 

i>i 

with the weights arranged in non-increasing order, w\ > W2 > . . . , and let us denote by 
Q — (£1 ■ (,i>)i t i>>i the matrix of scalar products of the points in the support of G. Let (cr l )i>i 
again be an i.i.d. sample from this measure and let Riy = cr 1 ■ cr 1 be the scalar product in 
£2 of cr 1 and cr 1 . For any n > 1 and a function / = /(cr 1 , . . . , er n ) of n configurations we will 
denote its average with respect to G®°° by 

(/>= E w h ---w ln f^ h ,...^ ln ). (1.4) 

«l,.-,in>l 

We will denote by E the expectation with respect to the randomness of G. 

In the mixed p-spin models, one expects the asymptotic Gibbs measures (jl.3p to be 
described precisely by the Parisi ultrametric ansatz (see [21], [II])- So far, most of the progress 
in the direction of proving structural results about G was based on the idea that certain 
information about its geometry can be recovered from the asymptotic stability properties of 
the Gibbs measure Gn under small perturbations of the Hamiltonian Hn{ct). In particular, 
two such stability properties in the setting of mixed p-spin model are well known - the 
Aizenman-Contucci stochastic stability [1] and the Ghirlanda-Guerra identities [13]. They 
can be written down in terms of the asymptotic Gibbs measure G as follows. 

1. (Ghirlanda-Guerra identities) Random measure G is said to satisfy the Ghirlanda-Guerra 
identities [13] if for any n > 2, any bounded measurable function / of the overlaps (Ri : i')i t i>< n 
and any integer p > 1 we have 

1 1 n 

IV ft 

1=2 

These constraints on the distribution of R look very mysterious but they arise from a very 
natural and very general principle of the concentration of the Hamiltonian Hn{ct) (see |13j). 
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2. (Aizenman-Contucci stochastic stability) Given integer p > 1, let (g p (£z))/>i be a Gaussian 
sequence conditionally on G indexed by the points (£z)z>i with covariance 

Cov(<7 p (6),<7 P (60) =(6 (1-6) 
Given tel, consider a new sequence of weights 

Wie tgpfo) 

wj = =— : -— T 1.7 

defined by a random change of density proportional to e t9p ^ l \ Let (wf) be the weights (it;*) 
arranged in the non- increasing order and let n : N — >■ N be the permutation keeping track of 
where each index came from, wf = w^ny Let us define by 

G* = <^(o and Q" = few ■ ^yoWi (1-8) 

the probability measure G after the change of density proportional to e t9p ^ and the matrix Q 
rearranged according to the reordering of weights. Measure G is said to satisfy the Aizenman- 
Contucci stochastic stability [1] if for any p > 1 and i el, 

(«)i>i, Q v ) = {{vn)t>i,Q) (1-9) 

where equality in distribution is in the sense of finite dimensional distributions of these arrays. 
This property represents the invariance of the distribution of measure G up to orthogonal 
transformations under the random changes of density (11. 7p and it arises from the continuity 
of the Gibbs measure Gn under small changes in the inverse temperature parameters (3 P in 
the Hamiltonian H N (cr) (see pQ). 

□ 

Originally, the proofs of the Ghirlanda-Guerra identities (jl.5p in |X 3j and the Aizenman- 
Contucci stochastic stability (11. 9p in [1] , [8] obtained these results for each p > 1 on average 
of the parameter ft p over any non-trivial interval. In a closely related formulation, one can 
always perturb the parameters ((3 p ) slightly (for example, one can find (Pn,p) such that 
\Pn, p -P p \< 2- p N- l ' l& ) so that the sequence of Gibbs measures Gn corresponding to these 
slightly perturbed parameters satisfies the above properties in the limit (see (33], [M])- Both 
of these formulations hold for mixed p-spin models with arbitrary subset of p-spin terms 
(II. ip present in the model, i.e. for which (3 P ^ 0. However, if the model contains terms for 
p = 1 and all even p > 2 then it was proved in [18] that the Ghirlanda-Guerra identities (II. 5p 
hold in the strong sense without perturbation of parameters for all integer p > 1. The proof 
was based on the validity of the Parisi formula for the free energy proved by M. Talagrand in 
[3Tj following the discovery of the replica symmetry breaking bound by F. Guerra in [12], on 
the differentiability properties of the Parisi formula (see [32], [17]) and on the positivity of 
the overlap (see [S])- A similar strong version of the Aizenman-Contucci stochastic stability 
was proved in [B]. 
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As we mentioned above, the importance of the stability properties (11. 5p and (11. 9p in 
the SK model comes from their many applications (see e.g. [I], [5], [6], [15], [16], [18], [2T] , 
[22] . [23] . [25] . [28] . [33] . [30] . [34]). Of course, the ultimate goal would be to prove the Parisi 
ultrametricity conjecture which states that the support of the asymptotic Gibbs measure G 
must be ultrametric in £2 with probability one, which would allow us to identify G with the 
Ruelle Probability Cascades in [27J. At the moment, the results that come closest to proving 
this conjecture are based either on the Aizenman-Contucci stochastic stability or on the 
Ghirlanda-Guerra identities. First such result was proved by L.-P. Arguin and M. Aizenman 
in [5] using the Aizenman-Contucci stochastic stability under a technical assumption that the 
scalar products £j • & of points in the support of measure G take finitely many non-random 
values. Following their work, a similar result was proved by the author in [15] (see also [22] 
for a recent elementary proof) and by M. Talagrand in [33] using the Ghirlanda-Guerra 
identities instead. Some modest progress toward the general case was made in [23] but the 
conjecture still remains open. Once this conjecture is proved, the Parisi formula for the free 
energy will easily follow using the Aizenman-Sims-Starr scheme developed in [2] which would 
naturally complete the mathematical justification of the Parisi ansatz in the SK model. It is 
worth mentioning that stability properties of the Gibbs measure under small perturbations 
of the Hamiltonian play very important role in other spin glass models as well (see e.g. [20]). 

In the main result of this paper we will show that one can combine the Ghirlanda- 
Guerra identities (II ,5p and the Aizenman-Contucci stochastic stability (II. 9p into a joint 
stability property as follows. It is known (Theorem 2 in [15]) that if the measure G satisfies 
the Ghirlanda-Guerra identities and if q* is the supremum of the support of the distribution 
of the overlap R12 under EG® 2 then with probability one G is concentrated on the sphere 
of radius y/q*. Let 

b p = ( q y-E(Rl 2 ). (1.10) 

Then the following holds. 

Theorem 1 A random measure G satisfies the Ghirlanda-Guerra identities U.5\) and the 
Aizenman-Contucci stochastic stability U.!J\) if and only if it is concentrated on the sphere 
of constant radius y/q* with probability one and for any p > 1 and t G M, 

(KW (9 P (^(i))-b P t) l > v Q n ) =(M*>r> (i-ii) 

where equality in distribution is in the sense of finite dimensional distributions. 

One can see that the Ghirlanda-Guerra identities are now replaced by the statement that 
the Gaussian field (fi'p(^)) after permutation ti corresponding to the reordering of weights 
in (II. 7p will only differ by a constant shift b p t in distribution. In the language of competing 
particle systems ([28] and [5]), (11.1 ip means that the past increments of the dynamics after 
re-centering have the same law as the future or forward increments. The stability property 
(11.1 ip is well-known for the ultrametric Ruelle Probability Cascades (see Theorem 4.2 in [7] 
or Theorem 15.2.1 in [31]) and, of course, the big question is whether it holds only for these 
measures and whether (11. lip implies ultrametricity. 
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The Ghirlanda-Guerra identities do not require the random measure G to be discrete 
and, in fact, the Aizenman-Contucci stochastic stability can be formulated not only for 
discrete measures as well. We will mention this more general formulation in the next section. 
However, we prefer to state our main result in the setting of discrete measures since it allows 
for a particularly attractive formulation (11. lip in the spirit of competing particle systems, 
as in [28] and [5]. Moreover, from the point of view of studying structural properties of such 
measures one can without loss of generality start with discrete measures since it is easy to 
show that sampling an i.i.d. sequence of points from the original measure and assigning them 
new independent weights from the Poisson-Dirichlet distribution creates a discrete measure 
which still satisfies both properties. On the other hand, almost any geometric property of 
the original measure will be encoded into a countable i.i.d. sample and, therefore, this new 
discrete measure. 

The unified stability property (II. lip inspired a new representation of the Ghirlanda- 
Guerra identities in [23] which yielded some interesting applications. For another recent 
stability property that reproduces the Ghirlanda-Guerra identities on average see [9] . 

Acknowledgment. The author would like to thank the referees for making many important 
suggestions that helped improve the paper. 

2 Proof of Theorem QJ 

Let (p l )i>i be an i.i.d. sequence from measure G n defined in (jl.8p and denote by Si t y = p l ■ p l 
the overlap of p l and p l . Analogously to (ll.4p . for any n > 1 and a function f = f(p l , . . . , p n ) 
of n configurations we will denote its average with respect to (G 71 )® 00 by 

(/>*= E W l--- W lf(^ihh---^Aln))- (2-1) 

h,...,l n >l 

We now will denote by E the expectation with respect to the randomness of G and the Gaus- 
sian sequence (g p ). Let us first make a simple observation that equality of finite dimensional 
distributions in ( II .9p and ( II. lip implies equality of averages with respect to the random 
measures in the following sense. 

Lemma 1 If U.ll\) holds then for any k > 1, any bounded measurable function f of the 
overlaps on k replicas and any integers ni, . . . , n& > 0, 

E(nfeo°') - Mr/covw)) = ^(u9p(^ i f(( R i,i'h'<k)) (2.2) 

l<k 77 l<k 

Under U.9]) . this holds with all ni = 0. 

Remark. One can consider ( 12. 2 p with all ni = as the definition of the Aizenman-Contucci 
stochastic stability for non-atomic measures in which case is the Gaussian field with 

covariance (11.61) . Moreover, in this case ( 12. 2 p should be considered as the analogue of (II. lip . 
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Proof. This is obvious by separating the sum in (11. 4p and (12.11) into finitely many terms 
corresponding to the largest weights and the remaining small weights. For example, 



l<k ' j\ j, 'I l<k 

h,-,3k<N l<k 

where the remainder IZ^ consists of the terms with at least one index ji, ■ ■ ■ ,j k > N. The left 
hand side of (12. 2p can be similarly broken into two sums. The finite sums are equal because 
they involve only finitely many elements of the arrays (II. lip which are equal in distribution 
by assumption. Thus, we only need to show that 1Z n becomes small for large N. First taking 
expectation in the Gaussian random variables (g p (£i)) conditionally on (wi) and Q and using 
that 

E(n>(4)P| W>Q) <L(n u ...,n k ) 

l<k 

we get that 



\Kn\< L(n 1 ,...,n k )\\f\\oo'& E w ^ "' w ^ - Lk E E 

(ji,-,jk<N) c j>n 



which goes to zero as iV — > oo. The remainder TZJj for the left hand side of (12.21) is controlled 
by exactly the same bound because, by (II. lip , the corresponding terms in 7V N and 71 n are 
equal in distribution. 

□ 

The "if part of the Theorem [1] is easy since assuming (11.111) we only need to prove (11.51) 
and this follows from integration by parts of (12.21) with n\ — 1, n-i — . . . — n k — 0. In this 
case the right hand side is zero by averaging g p (cr l ) first and the left hand side is 



E((^(p 1 ) - 6pf)/(0Mv'<*)) = m {(J2 s h ~ b p~ kS ik + i)f((Si,i'h'<*) 

i=i 

k 

= tE^R^i -K- kRi k+l )f({R hV \ v < k ) 
i=i 

k 

= tE((^Rli + mi2)-kRl k+1 )f((Ri,vh'<k) 



1=2 



where in the second line we used ( ll.9p part of (11.111) and Lemma [lj and in the third line we 
used (ll.lOp and the fact that = The fact that the last sum is zero is exactly ( II. 5p . 
To prove the "only if part we need the following key lemma. 

Lemma 2 // 1773]) and / fOj) hold then holds. 
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Proof. The proof is by induction on N = rii + . . .+rik- When N = 0, (I2.2p is the consequence 
of flUD by Lemma EJ Suppose (Q holds for all k > 1, all / and for all N < N . Clearly, 
we only need to prove the case of powers n\ + 1, n%, ■ ■ ■ , n^. Writing 

and using Gaussian integration by parts for g^cr 1 ) we can rewrite the right hand side of 
(12 .2p with rii + 1 instead of rii as 

Y^n^g^r . ..gpi^r- 1 ■ ..g P (<r k r k Rii /((^vv^)}- (2.3) 

l<k 

Again, writing 

(g^p 1 ) - = (^(P 1 ) " b p t){g p {p x ) ~ M) ni 

and using Gaussian integration by parts for g p {p 1 ) we can rewrite the left hand side of (12. 2p 
with rii + 1 instead of n\ as I + II where I is given by 

^n i E(fe(p 1 ) - b p t) ni . . . (g p (p l ) - b.t) 111 - 1 . . . (g p (p k ) - b p t) nk /((S w 0m'<*) ) (2-4) 
i<k 

and II is given by 

tE(n(^(pO-M) ni (E^-^-^ + i)/((V)M^))^ (2-5) 

l<k l<k ^ 

By induction hypothesis, (12 .4p is equal to (12. 3 p and (12. 5p is equal to 

tE(n^(^)"'(E i? M - b P - k Kk + i)f{^')^<k)). (2.6) 

l<k l<k 

Since (■) does not depend on the Gaussian sequence (g p (£i)) we can take expectation K g with 
respect to the randomness of this sequence conditionally on G first and notice that 

® 9 ng P (<T l T l = f , ((Ri,i>h><k) 

l<k 

for some function /' of the overlaps of k configurations cr 1 , . . . , er k . Therefore, (I2.6P equals 

k 

1 E ( (E R ii - b P~ kR lk + i) iff) ((%')«'<*) ) (2-7) 
i=i 

k 

= tE((j2Ki +W,2> - ^?, fe+ i)(//')((^ J /0^r< fe )) = 0, 

1=2 

where in the first equality we again used (11.101) and the fact that £j • £j = q* and the second 
equality is by the Ghirlanda-Guerra identities (II. 5p . This finishes the proof. 

□ 

The equality of joint moments (I2.2p proved in Lemma [2] implies the following. 
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Lemma 3 If U.5\) and U.9\) hold then 

{[g p {p l ) - b p t)^ v (<W>i) = (Ri,i'h'>i) ( 2 -8) 

where equality in distribution is in the sense of finite dimensional distributions. 

Remark. Let us recall that the i.i.d. sequences (cr l ) and (p l ) are sampled from G and 
correspondingly and, therefore, the distributions of the right-hand side and left-hand side in 
(El are under EG® 00 and E(G ,W )®°°. 

Proof. By choosing / to be monomials, (I2.2p gives the equality of joint moments of the 
corresponding elements of the two arrays in (12. 8p . In our case the joint moments uniquely 
determine joint distributions, for example, by the main result in [26] which states that we 
only need to ensure the uniqueness of one dimensional marginals and the fact that the one 
dimensional marginals are either bounded or Gaussian. 

□ 

Proof of Theorem [H Finally, we will show that (12. 8 p implies (II. lip . The procedure is 
very similar to the one at the end of Theorem 4 in [T5] or a more general argument in 
Lemma 4 in [19]. First of all, by the well-known result of Talagrand (Section 1.2 in [30]). 
the Ghirlanda-Guerra identities imply that the weights (wi) must have a Poisson-Dirichlet 
distribution PD(m) where m is determined by 

E(I(R lt2 = q*)) = E J2 w i = 1 - m - 

i>i 

This means that if (ui) is a Poisson point process on (0, oo) with intensity measure x~ m ~ x dx 
then wi — Ui/ Y2j>i u r ^ n particular, all the weights are different with probability one. This 
point is not crucial but it makes for an easier argument. The reason why (12.81) implies (11.111) 
is because one can easily reconstruct the arrays in (II. lip from the arrays (12.81) using that 
(a 1 ) is an i.i.d. sample from according to weights (wi) and (p l ) is an i.i.d. sample from 
(Ctt(z)) according to weights (wf). The key observation here is that given arrays (12. 8p we 
know exactly when a 1 = cr l and p l = p l since this is equivalent to R\y = q* and Sijr = q* . 
Therefore, given N > 1 and 

we can partition the set {l,...,iV} according to the equivalence relation I ~ /' defined 
by R[j> = q*, let the sequence of weights (wf)i>o be the proportions of the sets in this 
partition arranged in non-increasing order and extended by zeros and, given any integer 
j in the element of the partition corresponding to the weight , define ^ = aK We let 
Q N — (£/ V ' '£v)i,v>X- The elements of (^/ v ) and Q N with indices corresponding to zero weights 
can be set to some fixed values, and we break ties between wf by any pre-determined 
rule. Similarly, given 

((g P {p l ) ~ K^ikw ( S i,i')i,i'<n) 

we can construct sequences (wf), (£f) and Q N = (£f ■ £f). Equation (12. 8p implies that for 
any fixed k > 1, 

(WW (9 P an - v)K fc > (©w<*) =(wv, (9 P m\< h , Ki'h'<*)- 
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It remains to observe that the right hand side converges 



(M,<, (d50w<*) "^(Wk, Ww) (2-9) 



almost surely and, similarly, the left hand side converges a.s. to the corresponding array from 
the left hand side of (11. lip . To prove (12. 9L we notice that by construction 



is the empirical measure based on the sample a 1 , . . . , cr N from the measure G = ^ z>1 wrf^. 
By the strong law of large number for empirical measures (e.g. Theorem 11.4.1 in [H]), the 
laws G N — > G almost surely and since the Poisson-Dirichlet weights (wi) are all different a.s., 
the largest k weights must converge (wf)i<k —> (wi)i<k almost surely and for large enough 
N we must have (^/ v )z<fc = ((,i)i<k and, thus, (12. 9p holds. 
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